Two-dimensional dipolar Bose gas with the roton-maxon excitation spectrum 
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We discuss fluctuations in a dilute two-dimensional Bose-condensed dipolar gas, which has a roton- 
maxon character of the excitation spectrum. We calculate the density-density correlation function, 
fluctuation corrections to the chemical potential, compressibility, and the normal (superfluid) frac- 
tion. It is shown that the presence of the roton strongly enhances fluctuations of the density, and 
we establish the validity criterion of the Bogoliubov approach. At T = the condensate depletion 
becomes significant if the roton minimum is sufficiently close to zero. At finite temperatures ex- 
ceeding the roton energy, the effect of thermal fluctuations is stronger and it may lead to a large 
normal fraction of the gas and compressibility. 
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In the last decade, ultracold gases of dipolar particles, 
which include atoms with a large magnetic moment and 
polar molecules, attracted a great deal of interest [H-Il]- 
Being electrically or magnetically polarized such parti- 
cles interact with each other via long-range anisotropic 
dipole-dipole forces, which drastically changes the na- 
ture of quantum degenerate regimes. Experiments with 
chromium atoms (magnetic moment Q^jlb) 0]j together 
with theoretical studies d, |1] , have revealed the depen- 
dence of the shape and stability diagram of trapped dipo- 
lar Bose-Einstein condensates on the trapping geometry 
and interaction strength. They initiated spinor physics 
with quantum dipoles 0, Qi and now also dysprosium 
[o| and erbium [10] atoms (magnetic moments lO/is and 
7/LtB, respectively) have entered the game. Recently, fas- 
cinating prospects for the observation of novel quantum 
phases have been opened by the creation of ultracold 
clouds of polar molecules and cooling them to almost 
quantum degeneracy [1, [ll[. In this case, the dipole- 
dipole forces can be orders of magnitude larger, and one 
has a possibility to manipulate the molecules making use 
of their rotational degrees of freedom. 

For dipolar Bose-condensed gases, one of the key issues 
was related to the presence of the roton-maxon charac- 
ter of the excitation spectrum and to the possibility of 
obtaining supersolid states in which the condensate wave- 
function is a superposition of a uniform background and a 
lattice structure. The roton-maxon structure of the spec- 
trum has been first predicted for strongly pancaked Bose- 
Einstein condensates [12] ■ However, the idea to obtain a 
supersolid state when the roton touches zero and uniform 
BEG becomes unstable, did not succeed because of the 
collapse of the system [l^ Since that time, several 
proposals have been made for the creation of supersolid 
states with bosons . They rely on the potential of in- 
teratomic interaction which is flat at short distances and 
decays at large separations, and the results then indicate 
the presence of dense supersolid clusters [T6| . 

This activity brought in analogies with liquid helium, 



where the studies of the roton-maxon spectrum and the 
attempts to observe experimentally the supersolid state 
have spanned for decades [O, ESi] after the early theo- 
retical prediction [l^. However, the most credible claim 
for the observation of the supersolid [13] is now with- 
drawn It is worth mentioning that the old idea of 
obtaining a stable density-modulated state (supersolid) 
in superfluid helium moving with a supercritical velocity 
p^ l is now being discussed in a general context of super- 
fluidity in Bose gases flowing with velocities larger than 
the Landau critical velocity |23| . 

The roton-maxon character of the excitation spectrum 
was also attracting large attention by itself. In relation 
to liquid helium, it has been discussed how the position 
of the roton minimum influences the phenomenon of su- 
perfluidity fis!]. In the context of dipolar bosons in two 
dimensions, numerical calculations of zero temperature 
phase diagram (23 - l26| found that the reduction of the 
condensed fraction with an increase in the density can 
be attributed to the appearance of the roton minimum. 
Finite temperature Monte Carlo calculations [13] have 
revealed that the rotonization of the spectrum can de- 
crease the Kosterlitz-Thouless superfluid transition tem- 
perature. Although the calculations ^3 - l27| were focused 
on fairly high densities, they raised the question of ap- 
plicability of the Bogoliubov approach for dipolar bosons 
|26{ . At the same time, there was an extended activity 
on the static and dynamical properties of dilute trapped 
dipolar Bose-Einstein condensates on the basis of the Bo- 
goliubov approach [28l - l3l| . 

It is therefore instructive to identify the validity cri- 
terion of the Bogoliubov approach for Bose-condensed 
dipolar gases with the roton-maxon excitation spectrum, 
and this is the subject of the present paper. We show 
that at zero temperature the density fluctuations orig- 
inating from the presence of the roton minimum, lead 
to a signiflcant depletion of the condensate and modify 
thermodynamic quantities if the roton minimum is suf- 
flciently close to zero. At flnite temperatures exceeding 
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the roton energy, thermal density fluctuations may have 
a much stronger influence, leading to a large increase of 
the normal fraction and compressibility. 

We consider a dilute Bose-condensed gas of dipolar 
bosons (tightly) confined in one direction (z) to zero 
point oscillations and assume that in the x,y plane the 
translational motion is free (see Fig.l). The dipole mo- 
ments are oriented perpendicularly to the x,y plane, 
which for electric dipoles (polar molecules) can be done 
applying an electric field, and for magnetic atoms by us- 
ing a magnetic field. In this quasi-2D geometry, at large 
interparticle separations r the interaction potential is 



(1) 



with d being the dipole moment, m the particle mass, 
and = mcP/h^ the characteristic dipole-dipole dis- 
tance. The short-range part of the potential is assumed 
to be such that there is a roton-maxon excitation spec- 
trum. In the ultracold limit where the particle momenta 




FIG. 1. Dipolar Bose-Einstein condensate tightly confined in 
one direction. 

satisfy the inequality fcr* 1, the off-shell scattering am- 
plitude defined as f{k,k') = J exp{—ik'f')V{r)^^{r)cPr 
(ipki^ is the wavefunction of the relative motion with 
momentum fc), is given by (see [32| and refs. therein): 
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where we omit higher order terms in fc. The second term 
in the right hand side represents the so-called anomalous 
contribution coming from distances of the order of the de 
Broglie wavelength of particles fsS*] . It takes into account 
all partial waves and is obtained using a perturbative ap- 
proach in V{r). The first term in the right hand side of 
Eq.(I2]) describes the short-range contribution. It is ob- 
tained by putting k' = and proceeding along the lines of 
the 2D scattering theory [s^]- The parameter k depends 
on the behavior of V{r) at short distances. In the quasi- 
2D geometry it also depends on the confinement length 
in the z-direction, Iq = yJh/mujQ, where Wg is the con- 
finement frequency. One then can express k through the 



3D coupling constant ga^). If the 3D s-wave scattering 
length, a^p = mg^D/ATrh^ Iq, then k is exponentially 
small [3J, |35| and we may omit the k-dependence under 
logarithm in Eq.(I2|), as well as i7r/2 in the denominator of 
the first term. This gives /(fc, fc') = g{l-'C\k~k'\), where 
the 2D short-range coupling constant is g ^ gso/V^lo 
and C = 2TTh^r^/mg = 2ttcP /g. Employing this result in 
the secondly quantized Hamiltonian , we obtain 
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where S is the surface area, = ?i^fc^/2m, and a]^, 
Ok are the creation and annihilation operators of parti- 
cles. At zero temperature there is a true Bose-Einstein 
condensate in 2D, and we may use the standard Bogoli- 
ubov approach. Assuming the weakly interacting regime 
where 'mg/2'KTT? <C 1 and r* <C ^, with ^ = h/ ^Jrrvng 
being the healing length, we reduce the Hamiltonian 
(j3]) to a bilinear form, use the Bogoliubov transforma- 
tion at — Ukb^j^ — Vkb_^, and obtain the diagonal form 

H = Eq + X^fe^fe^l^fe ill terms of operators bt, of el- 
ementary excitations. The Bogoliubov functions Uk,Vk 
are expressed in a standard way: Uk,Vk ~ {^Sk/Eu ± 

Ek/ ek) /2, and the Bogoluibov excitation energy is 
given by Efc = ^ E^ + 2ngEk{l - Ck). 

To zero order the chemical potential is fi = ng. For 
small momenta the excitations are sound waves, Ek — 
^ ng/mk. The dependence of ek on fc remains monotonic 
with increasing fc if C < -\/8^/3 (see Fig. 2 ). For the 
constant C in the interval 
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the excitation spectrum has a roton-maxon structure. It 
is then convenient to represent Sk in the form: 
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where fc^ = 2C /^ and fcA = \/4/C^ — fc^. If the roton is 
close to zero, then fc^ is the position of the roton, and 



/S. — 1i^krk/^/2m — 2ngC\Jmng/h^ — C^^mng/h^)^ (6) 

is the height of the roton minimum (see Fig. 2). For C — £, 
the roton minimum touches zero, and at larger C the 
uniform Bose condensate becomes dynamically unstable. 

It should be noted that the coupling constant g can 
be tuned by using Feshbach resonances or by modifying 
the frequency of the tight confinement wq- Therefore, 
although the range of C given by Eq.Q is rather nar- 
row, it can be reached without serious difficulties. The 
condition C = 2nd^/g = £, is reduced to {mg/2iTfi^) = 
asD / V^Iq — 2nnrl . For dysprosium atoms we have the 
dipole-dipole distance ~ 200 A, and at 2D densities 

10^ cm~^ the roton-maxon spectrum is realized for the 
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3D scattering length a^o of several tens of angstroms at 
the frequency of the tight confinement of 10 kHz leading 
to the confinement length Iq about 1000 A. 




FIG. 2. Excitation energy Sk of the quasi-2D dipolar BEG 
as a function of momentum k for several values of kr- The 
solid curve (krS, = 1.84) shows a monotonic dependence Sk, 
the dotted curve (fc^f = 1.96) is Sk with the roton-maxon 
structure, and the dashed curve (fer^ = 2.08) corresponds to 
dynamically unstable BEG. 

The Bogoliubov approach assumes that the density 
and phase fluctuations are small. In the 2D case at T = 0, 
the presence of the roton does not significantly change the 
phase fluctuations and they remain small. However, the 
situation with the density fluctuations is different. Writ- 
ing the operator of the density fluctuations as (see (36j ) 

5h ~ \/nJ2-kiiuk — Vk) exp(jkr)6k + /i-c., we obtain for the 
density-density correlation function: 
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where — [exp(ej,/T) — 1]^^ are occupation numbers 
for the excitations. It is instructive to single out the roton 
contribution to the correlation function ([T]). Asuming 
that the roton is close to zero and the roton energy is 
A ^ ng, we have the cofficient C close to ^, and kr ~ 
Then, using Eqs.([5]) and (O, for the contribution of 
momenta near the roton minimum at T = we obtain: 

{6h{r)5n{0)) r ^mg f 2ng 



In 



] M2r/0; A « ng, (8) 



where Jq is the Bessel function. 

We thus see that the density fluctuations grow log- 
arithmically when the roton minimum is approaching 
zero and they can become strong for very small A. 
In this case they lead to a significant depletion of the 
condensate. The non-condensed density of particles is 
n' — J vf.(P k / s-nd the integral over dk is logarith- 
mically divergent at large momenta because of the dipo- 
lar contribution to the interaction strength, —gCk. How- 
ever, this form of the dipole-dipole contribution is valid 
only for k <C l/r*. We thus may put a high momentum 
cut-off 1/r*, which leads to (see Fig. 3): 
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In the absence of the dipole-dipole interaction (r* = 
and kr — 0) we recover the usual result for the 
2D BEC with short-range interparticle repulsion, n' — 
n{mg / AttTl^) . For A < we have (2 - krS) ~ (fcAO^/4 
and Eq.Q transforms to 
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where C = \j li^fi' j e^mg . 

As we see from Eqs. ([8]) and ([T0| . for the roton min- 
imum close to zero a small condensate depletion and 
small fluctuations of the density require the inequality 
{mg I 'Kh^)\vL{2ng 1 1\) <C 1. It differs only by a logarith- 
mic factor ln(2n(7/A) from the small parameter of the 
theory, {mg/2'KTi^) <^ 1, in the absence of the roton. 
The same logarithmic factor appears in the fluctuation 
correction to the chemical potential and in the one-body 
density matrix gi{r) — (5'^^(r)^(0)), where vE'(r) is the 
field operator. Assuming that the roton minimum is close 
to zero and taking into account only the contribution of 
momenta near this minimum we have for A <^ ng: 



gi{r) = no 
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The correction to the chemical potential due to quantum 
fluctuations is given by: 



(5/i 



2m g 



A < ng. 



(12) 
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FIG. 3. Non-condensed fraction as a function of kr^ for 
mg/Aivh^ = 0.01 (^/r, — 100/fcr^). A similar increase of 
the non-condensed fraction with decreasing the roton energy 
A has been found in numerical calculations of Ref. [2811. 



However, the situation changes in the calculation of 
the compressibility. At T = the inverse compressibility 
is equal to n^d^/dn. Then, using Egs. p^ and ([B]), for 
the roton minimum close to zero we obtain at A <C ng: 



dn ^ 



2mg ( 2ng\ mg ( 2ng\ 



(13) 



where g is the mean field contribution, and the second 
and third terms originate from quantum fiuctuations. 
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Small deviations of the compressibility from the mean 
field result require the inequality 



quantity reads (c.f. [36 



< 1. 



(14) 



We thus conclude that at T = the validity of the Bo- 
goliubov approach is guaranteed by the presence of the 
small parameter (jl4p . For the dysprosium example given 
after Eq.® we have ng about 5 nK, and the criterion 
([14]) is satisfied for the roton energy above 2 nK. 

In 2D at finite temperatures, long- wave fluctuations 
of the phase destroy the condensate [37l - [39j . There is 
the so-called quasicondensate, or condensate with fluc- 
tuating phase. In this state fluctuations of the density 
are suppressed but the phase still fluctuates. The tran- 
sition from a non-condensed state to quasiBEC is of the 
Kosterlitz-Thouless type and it occurs through the for- 
mation of bound vortex-antivortex pairs ^40] . Somewhat 
below the Kosterlitz-Thouless transition temperature the 
vortices are no longer important, and in the weakly in- 
teracting regime that we consider the phase coherence 
length is exponentially large. Thermodynamic proper- 
ties, excitations, and correlation properties on a distance 
scale smaller than are the same as in the case of a 
true BEC. Moreover, for realistic parameters of quantum 
gases, exceeds the size of the system [4lj, so that one 
can employ the ordinary BEC theory. 

Irrespective of the relation between and the size of 
the system, one may act in terms of the density and phase 
variables (hydrodynamic approach). We now show that 
the rotonization of the spectrum can strongly increase 
thermal fluctuations of the density and destroy the Bose- 
condensed state even at very low T. Using equation ([7|) 
we calculate the density-density correlation function. As- 
suming that the roton energy A is very small (at least 
A ^ T), the main contribution to the integral in Eq.© 
comes from momenta near kr, and we obtain: 



((5n(r)(5(0)) _ Amg T 

^2 - 
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where it is also assumed that k^r ^ 1. Comparing this 
result with Eq. ^ we see that instead of the logarithmic 
factor we have 2ttT/A > 1. 

The same factor appears in the correction to the chem- 
ical potential due to thermal fluctuations: 



(5/i 
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We now calculate the density of the normal component 
in the presence of the roton. In 2D the expression for this 



riT = 



h^k"^ dNk d^k 



2m dek (27r)2 ' 
If the roton minimum is close to zero and A ^ T, then 
the momenta near the roton minimum are the most im- 
portant, and the integration over dk yields: 
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The employed approach requires the condition ut <C n 
because we used the spectrum of excitations obtained by 
the Bogoliubov method. Again, at temperatures T > A 
we should have the inequality {2mg/h'^)T/A <^ 1. 

A different small parameter appears in the calcula- 
tion of the compressibility. The inverse isothermal com- 
pressibility is proportional to {dP/dn)T, where the pres- 
sure is _P = —{dF/dS)T, with the free energy given by 
F = Eo +TY,i,H^ - exp{-ek/T)]. For the roton mini- 
mum close to zero and {T, ng} 3> A, we obtain: 



dP 
dn 



ng 



mg f 2ng\^ T 
n'' \A ) A 
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where we omitted less important finite-temperature con- 
tributions and the zero temperature contribution propor- 
tional to the small parameter ([T4| . Eq. ([T8|) shows that at 
T ^ A the Bogoliubov approach requires the inequality 



(19) 



whereas for T < A it is sufficient to have criterion ([T] 

For certain quantities the Bogoliubov approach may 
give good results at T > A if {nig/n^)T/A < 1, and at 
T < A in the presence of the ordinary small parameter 
nng/2'Kh^ amplified by a logarithmic factor ln(2ng/A) for 
ng 3> A. However, the validity of this approach is guar- 
anteed only if the inequalities ([T^ and (|14p are satisfied. 
For T ^ A the compressibility following from Eq. p^ 
and the normal fraction given by equation (jl7p can be- 
come significant, by far exceeding similar quantities in 
an ordinary 2D Bose gas with short-range interactions at 
the same temperature, coupling constant g, and density. 

In conclusion, we have shown that the roton-maxon 
structure of the excitation spectrum, which can be 
achieved in dipolar Bose-condensed gases in the 2D ge- 
ometry, strongly enhances fluctuations of the density. 
We obtained the validity criterion of the Bogoliubov ap- 
proach and found that at finite temperatures in the di- 
lute regime where nr^ <C 1, thermal fiuctuations may 
significantly increase the compressibility and reduce the 
superfluid fraction even at temperatures well below the 
Kosterlitz-Thouless transition temperature. 
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